A parameterisation of single and multiple muons in 
the deep water or ice 



A Margiotta 

Dipartimento di Fisica dell'Universita and Sezione INFN di Bologna, 
viale C. Berti-Pichat, 6/2, 40127 Bologna, Italy 

E-mail: Annarita.Margiotta@bo . infn. it 

Abstract. A new parameterisation of atmospheric muons deep underwater (or ice) is 
presented. It takes into account the simultaneous arrival of muons in bundle giving the 
multiplicity of the events and the muon energy spectrum as a function of their lateral distribution 
in a shower. 



1. Introduction 

Primary cosmic rays interacting with atmospheric nuclei produce electromagnetic and hadronic 
showers that, during their development through the atmosphere, originate a large number of 
neutral and charged particles. Although neutrino telescopes are located under large depths of 
water or ice, a great number of high energy atmospheric muons can reach the sensitive volume 
of the detectors [Q. Indeed, they are the most abundant signal in a neutrino telescope and can 
be used to calibrate the detector and to check its expected response to the passage of charged 
particles. On the other hand, they can represent a dangerous background source for neutrino 
telescopes if track reconstruction programs, which use Cherenkov light emitted by downgoing 
muons, yield an upward going muon track [2J. 

A full Monte Carlo simulation, starting from the simulation of atmospheric showers, can 
accurately reproduce the main features of muons reaching a neutrino telescope, but requires 
a large amount of CPU time. Some fast parameterisations of the atmospheric muon flux 
underwater are available in the literature [3J @J [3] but all of them do not take into account the 
simultaneous arrival of muons in bundles, which are particularly dangerous for misreconstructed 
tracks. 

In this work, parametric formulae are presented to evaluate the flux of atmospheric muons, 
taking into account the muon multiplicity and the muon energy spectrum in a bundle, as a 
function of the distance from the shower axis. A detailed description can be found in jSj. 



2. The Monte Carlo simulation 

The parameterisation of the multiple muon flux and energy spectra presented here relies on a 
full Monte Carlo simulation of primary Cosmic Ray (CR) interaction and shower propagation 
in the atmosphere (HEMAS code [7j). The adopted primary CR flux and composition model is 
described in (Hj. The muons from the decay of secondary mesons reaching the sea level are then 
propagated to seven depths, from 2.0 kmw.e. down to 5.0 kmw.e., in steps of 0.5 kmw.e., using 



the MUSIC code jSj. The energy spectrum of muons depends on the vertical depth h, on the 
zenith angle 9, on the muon multiplicity in the shower m and on the distance of the muon from 
the shower axis R. Starting from the distributions of muon multiplicity, energy and distance of 
muons from the shower axis, a number of fits was performed in order to obtain the values of the 
constants to be used in the analytic description of the muon flux and energy. 



3. The parameterisation 

• The multiplicity distribution is reproduced, as a function of water depth and zenith angle, 
following an analytic description used by the Frejus collaboration [TUj : 
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A ~ is compatible with all our fits. 9 constants are necessary to describe eq. \I\m terms of h 
and i?. From (^Q), the depth-intensity relation for vertical muons (i? = 0°) can be calculated. 
A comparison with other parameterisations is shown in fig. 1. In fig. 2 comparisons of the 
zenith distribution evaluated with eq. ^and some experimental data are shown. 
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Figure 1. Muon vertical intensity versus 
depth. Comparisons are shown among 
this parameterisation for single muon events 
(dashed line) and for the total muon flux 
(full line), other parameterisations: Okada [3] 
(dotted-dashed line), Bugaev et al. Hj (dotted 
line) and the full Monte Carlo simulation 
(points). 



Figure 2. Muon flux versus zenith angle. 
Dot-dashed line: muon flux evaluated at 
h = 1.64 km w.e. compared with AMANDA- 
II measurements (triangles). Full line: 
muon flux evaluated at h = 3.8 km w.e. 
compared with the angular distribution (1- 
sigma error) measured by the NESTOR 
Collaboration ^2] (dashed lines). 



The energy spectrum of single muons is described by: 
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(|13j). where 7 is the spectral index of the primary beam and e = a/j3 is defined from 
~{ ~Ex^ ) = a + (3E^. In this work, 7 and e were used as fit parameters depending on h 
and i?; (5 was fixed to a value of 0.420 (kmuu.e.)" 1 . 6 constants are necessary in eq. |2 
The muon lateral spread in muon bundles is parameterised as 

^ = C R (3) 
dR (R + R Q ) a {1 



R and a can be expressed as functions of h, m and $ using 9 constants. The description 
of the muon lateral spread is the preliminary step to evaluate the muon energy distribution 
in a bundle. 

• The energy spectrum of muons arriving in bundles has the same general form as for single 
muons, but each parameter shows a dependence on multiplicity and distance from the axis 
and not only on depth and zenith angle. 15 constants are necessary in this case. 

A comparison between the average energy of single and double muons evaluated with 
this parameterisation and the results of MACRO measurements, performed with a Transition 
Radiation Detector |14j under different depths of rock, is presented in Tab. 1. A correction has 
been applied to take into account the differences in energy losses of muons in standard rock and 
water. 
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Table 1. Average energy of single (column 4) and double (column 6) muons measured by the MACRO 
collaboration |14j for different values of the slant depth (column 1) and of the average muon zenith 
angle direction (column 2). The predictions of this work are given in columns 5 and 7 for m = 1 and 
r7i = 2, respectively. The calculated values have been corrected to take into account the muon energy 
loss differences in water and rock. For double muons, the computed average distance of muons from the 
shower axis, reported in column 3, has been assumed. 



4. Conclusions 

A new parameterisation of the atmospheric muon flux under water, considering the simultaneous 
arrival of muons in bundles, is presented. Energy and zenith distributions obtained with these 
formulae are in good agreement with the predictions from other parameterisations, with the full 
Monte Carlo and with experimental measurements. This parameterisation can be used as a fast 
generator to simulate atmospheric muon fluxes in underwater /ice neutrino telescopes. 
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